Nonequilibrium Josephson effect in short-arm diffusive SNS interferometers. 
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We study non-equilibrium Josephson effect and phase-dependent conductance in three-terminal 
diffusive interferometers with short arms. We consider strong proximity effect and investigate an 
interplay of dissipative and Josephson currents co-existing within the same proximity region. In 
junctions with transparent interfaces, the suppression of the Josephson current appears at rather 
large voltage, eV ~ A, and the current vanishes at eV > A. Josephson current inversion becomes 
possible in junctions with resistive interfaces, where it occurs within a finite interval of the applied 
voltage. Due to the presence of considerably large and phase-dependent injection current, the critical 
current measured in a current biased junction does not coincide with the maximum Josephson 
current, and remains finite when the true Josephson current is suppressed. The voltage dependence 
of the conductance shows two pronounced peaks, at the bulk gap energy, and at the proximity gap 
energy; the phase oscillation of the conductance exhibits qualitatively different form at small voltage 
eV < A, and at large voltage eV > A. 

PACS numbers: 74.25.Fy, 74.45. -fc, 73.23.-b 



I. INTRODUCTION 

Multiterminal superconductor - normal metal - super- 
conductor (SNS) junctions are interesting devices where 
an interplay between the dissipative normal electron cur- 
rent and non-dissipative Josephson current can be stud- 
ied. The simplest device of this type consists of two su- 
perconducting reservoirs and one normal reservoir con- 
nected by a small normally conducting T-shape bridge, 
see Fig.^a). A mesoscopic size of the bridge is essential 
to keep the coherence of the current transport over the 
whole device. During the last decade, a large amount of 
interesting experiments have been done with such kind of 
devices (for the review see Ref. 1 and references therein, 
further references can be found in Ref. 0) . 

Non-equilibrium state in multiterminal SNS junctions 
exhibits two closely related major phenomena: the inter- 
ferometer effect, which concerns the dependence of nor- 
mal conductance of the device on the phase difference 
between the superconducting reservoirs?^ and the Joseph- 
son transistor effect, which concerns the dependence of 
the Josephson current on the current injected from the 
normal reservoir^i^ The interferometer effect gives rise 
to a number of so called Josephson-like effects&i 

The interferometer effect has received the most of at- 
tentions; it has been extensively studied experimental- 
lyL^ and theoreticallyfiiSiS and presently this effect is 
rather well understood. The proximity of the supercon- 
ducting reservoirs leads to a modification of the density 
of states and transport properties of the normal bridge 
(proximity effect), which therefore become sensitive to 
the phase difference at the reservoirs, and exhibit oscil- 
lating behavior as the function of the phase difference. 

The Josephson transistor effect has been also observed 
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feet has been first predicted^ and then theoretically stud- 
ied in the ballistic junctionsiSiiiiiii^ Non-equilibrium 



FIG. 1: Sketch of 3-terminal SNS interferometer: T-shape 
diffusive metallic wire with short arms (shaded region) is con- 
nected to superconducting (S) and normal (N) reservoirs (a), 
flux bias setup for measuring current-phase relation (b), cur- 
rent bias setup for measuring critical current (c). 

population of the Andreev statesi^ induced by the cur- 
rent injection leads to a full-scale variation of the Joseph- 
son current with the applied voltage and to the inversion 
of the current direction. Similar transistor effect has been 
also investigated in diffusive j unctions piiiiSiiSiSfi 

The most of experiments with multiterminal junctions 
have been done using diffusive metallic bridges, whose 
length was large compared to the superconducting co- 
herence length ^0- III such long junctions, the proxim- 
ity effect is suppressed, which results in small ampli- 
tude of the conductance oscillation (typically few per- 
cents of the full conductance value) , and small magnitude 
of the Josephson current. In theoretical analysis of the 
interferometer effect in diffusive junctions, the presence 
of the Josephson current is usually ignored. Similarly, 
the theoretical studies of the nonequilibrium Josephson 
g-g:g(,^i7ji8ji9 g^j.g j-estricted to the regime of "weak" prox- 
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imity effect jSiSiiSS when the induced gap in the normal 
bridge is much smaller than those in the superconduct- 
ing reservoirs. Such a regime is relevant for long diffu- 
sive junctions and for junctions with high-resistance NS 
interfaces. In addition, the non-equilibrium Josephson 
current was calculated for specific four-terminal circuit 
geometries where the dissipative normal current and the 
Josephson current were spatially separated. 

Meanwhile, it is conceptually interesting to investigate 
the problem of interplay of dissipative and Josephson cur- 
rents flowing in the same diffusive lead under a strong 
proximity effect. In this paper, we address this prob- 
lem by studying the non-equilibrium Josephson effect in 
3-terminal SNS interferometers with short arms, whose 
lengths are smaller than the coherence length ^o- 

Co-existence of dissipative and non-dissipative cur- 
rents in the same proximity region makes it difficult to 
identify the Josephson current component. In equilib- 
rium, the total current flowing through a proximity lead 
is the Josephson current, which is entirely determined 
by the supercurrent spectral density and the population 
numbers of the relevant states. The presence of the nor- 
mal injection current and related gradients of the distri- 
bution functions violates the local conservation of the su- 
percurrent component which varies along the lead, and 
hence its direct connection to the Josephson current is 
lost. Nevertheless, as we will show in the paper, a simple 
picture of the non-equilibrium Josephson effect as the 
result of non-equilibrium population of the states with 
the same current spectral density as in equilibrium can 
be justified for some particular cases (see also Ref. IT^ . 
However, a general situation seems to be more complex. 

In the short-arm SNS interferometers, the proximity 
effect is strong when the interfaces are transparent, but 
it also may be strong when the interfaces are highly resis- 
tive. The measure of the strength of the proximity effect 
is the magnitude of the induced energy gap, which is 
comparable with the superconducting energy gap in the 
reservoirs. Consequently, there are full-scale variations of 
the Josephson current and the normal conductance with 
the applied voltage and the phase difference. 

In short diffusive junctions with transparent interfaces, 
the Josephson current is solely carried by the Andreev 
states whose energies are smaller than the supercon- 
ducting energy gap, which is similar to the short bal- 
listic junctions. ^'^'^ Consequently, in these junctions, the 
Josephson current can be suppressed to zero but never be 
reversed. The Josephson current reversion becomes pos- 
sible in the junctions with resistive interfaces due to the 
negative contribution of the states with energies above 
the bulk energy gap. Moreover, in contrast to weak prox- 
imity regime^ the current reversion exhibits a fine struc- 
ture similar to the one theoretically discussed for ballistic 
junctions, and recently demonstrated in the experiment 
with long diffusive junctions^ 

The interplay of the injected and Josephson currents 
in the strong proximity regime is important for the inter- 
pretation of experiments with current biased junctions. 



Fig. ^c) . The magnitude of the non-equilibrium critical 
Josephson current measured with this setup is different, 
as we will show, from the one derived from the current- 
phase relation measured by an rf-SQUID, Fig.^b). This 
phenomenon has been earlier noticed in the absence of 
true Josephson current in the weak proximity regime^ 

The structure of the paper is the following. After 
introducing a basic formalism in Section II, we discuss 
the spectral functions in Section III, and the distribution 
functions in Section IV. Sections V and VI are devoted 
to a discussion of the non-equilibrium Josephson effect; 
the interferometer effect is considered in Section VII. 



II. BASIC EQUATIONS 

The junction we are going to investigate is sketched in 
Fig. ^a). It consists of two superconducting reservoirs 
and a normal reservoir connected by mesoscopic T-shape 
diffusive metallic bridge. Such a geometry can be realized 
in experiment, e.g., by using nanowire technology.^"* The 
superconducting reservoirs are assumed to have equal po- 
tentials, and the superconducting phase difference be- 
tween the reservoirs is The distance between the 
superconductors, 2d {—d < x < d), and the length of 
the injection lead, L {0 < y < L), are assumed to be 
small compared to the superconducting coherence length 
^0 — y^TiV/ A (T) is the diffusion coefficient), however, 
the relation between these lengths can be arbitrary. For 
simplicity, we assume the cross sections and normal con- 
ductivities of all wires to be equal, and the current from 
the voltage-biased normal reservoir to be injected in the 
middle of the SNS junction. 

Neglecting spatial variations of all quantities across the 
leads, we use one-dimensional static equations^ for the 
4x4 matrix Keldysh-Green function G in the normal 
leads, in which we neglect the inelastic collision term, 

[a.E, G] = ihVdJ, J = GdG, & = i, (1) 
^"(^o" ?^)' G'^^g^f-fY- (2) 

Here g^'"^ arc the retarded and advanced Green func- 
tions, f — f+ +azf- is the matrix distribution function, 
and d denotes spatial derivative. At the junction node, 
the matrix current J obeys the Kirchhoff's rule^S^ 

Jx=-0 = Jx=+o + Jy=+o- (3) 

The Keldysh component of the matrix current J de- 
termines the electric currents in the leads, 

I'OO I-OO 

j = f/ dETraJ'' = - dE I^{E), (4) 
4e Jo e Jo 

where a is the normal conductivity. The current spectral 
density I^{E) in Eq. 10} has three componcnts^^ 

/_ = (1/4) Tra, = + IJ+ - Iandf+, (5) 



3 



The first term in Eq. jSj) describes a dissipative current 
which provides usual Drude conductivity in the normal 
state. The second term gives conventional Josephson cur- 
rent in equilibrium, while the third term, the anomalous 
current, appears in non-equilibrium superconducting 
j mictions. Another diagonal component of J^, 

/+ = (1/4) Tr = D+dU + Isf- + landf-, (6) 

has the meaning of the net quasiparticle current (the sum 
of the electron and hole probability currents). Explicit 
equations for the spectral characteristics of the junction, 
D± , Is , and Ian , are conveniently written in terms of the 
following parametrization of the matrix g, 



g — u + V ^ a^u + exp{iaztp)iayV, 



1, (7) 



The function u is related to the quasiparticle density of 
states (DOS), normalized by its value in the normal state, 
N{E,x) = Reu^{E,x), while the function v is related 
to the spectral density of the condensate. The complex 
phase ip appears in the presence of a supercurrent. In 
these notations, the diffusion coefficients read 



D± = (l/4)Tr(l -u^u^T^)^*^) (8) 
=F \ v"\ ' cosh (2Im'0^ 



(1/2) 1 + lu^l^ - L.«P„.„i, /oT™„;,fl^ 



In the normal state, D± turn to unity. 

The spectral densities of the supercurrent and ano- 
malous current are given by equations. 



Is = (l/4)Trcr^ {v^dv'^-v'^dv'^) = - Im (i^^ai/') , (9a) 
Ian = (1/4) Trcr^-y^f)-^ = - \v"f sinh (2ImV'^) /2.(9b) 



In Eqs. © and ®, the relations (w, u)-^ = -{u,v)^* 
and ip'^ = ip^* are used, which follow from the general 
relationship g^ — —a^g^^crz-^^ 

Calculation of the electric current in Eq. involves 
the two steps: first one has to solve the Usadel equations 
for the Green functions g^'^, and then to solve the kinetic 
equations to find the distribution functions. 



III. SPECTRAL FUNCTIONS 

The Green function components of Eq. ^ represent 
the Usadel equations for the spectral functionS)2i 

2Ev ^ ihV [d{udv - vdu) - uvidtp f] , (10) 
v^dtp = I, (11) 

where the spatial constant I{E) is related to the super- 
current spectral density in Eq. Is — — Im/^. In 
terms of the spectral angle 9 related to the spectral func- 
tions as u = cosh0, V = sinh0, Eq. H1U|) takes the form. 



The two terms in r.h.s. of Eq. (|12(l are related to two dif- 
ferent depairing mechanisms, which provide spatial de- 
crease of 9 towards the middle of the junction. The first 
term is associated with the dephasing between the elec- 
tron and hole wave functions at finite energy E. The sec- 
ond term describes depairing caused by the time-reversal 
symmetry breaking due to supercurrent flow. 
The solution of Eqs. ((TT|l and ifT^ reads 



KE.x,) 



'iP{E,x,)^tPo+I 



--, Xi = x,y, (13) 
^ dz 

/o ^;%b:i)' ' 

R{E,9) = C+{4:E/ihV) cosh 9- (//sinh6l)^ (15) 

where C{E) is the integration constant, and V'o are 
the spectral functions at the junction node. The bound- 
ary conditions for Eqs. (|ll|l and H12|l are imposed by the 
conservation law for the matrix current in Eq. 



r=-0 = 9^0\x=+O + 9: 



(16) 



and similar equation for the function ip. At the normal 
electrode, the current I turns to zero, which means that 
7=0 along the injection lead, and therefore tp{y) = 
const. Hence, the derivative of ■0 is continuous in the 
horizontal lead, whereas 9{x) has a kink at a; = 0. In 
what follows, we consider symmetric junction, in which 
9{x) is even, and tp{^) is odd function. In this case, the 
phase tpi together with the anomalous current Ian, turns 
to zero at the junction node, and the kink in 6{x) is 
symmetric, 9'^ = d^Q\.j:=+fi = -dxQ\x=-ii- 

The boundary conditions at the NS interfaces depend 
on the interface resistance. Below we analyze the two 
different situations related to perfect and high-resistive 
interfaces, respectively. 



A. Transparent interfaces 

If the interface electric resistance Rt^s is much smaller 
than the normal resistance Rjsi of the horizontal lead, 
one can assume the spectral functions to be continuous 
at a; = id, namely, ^/'(±d) = ±0/2, and 0(±d) — 9s = 
Arctanh (A/i?); at the normal electrode 6* = 0. In this 
limit, the second term in the r.h.s. of Eq. H15() can be 
neglected due to large gradients of the spectral functions 
along the leads, which results in linear change of 9 along 
the injection lead, 9(E^ y) = 9o{l — y/L), and the bound- 
ary condition in Eq. (|16|l takes the form 



(17) 



d^e = {2E/iKD) sinh 9 + 1"^ cosh 61/ sinh^ 9. 



(12) 



The analytical expressions for the spectral functions in 
the horizontal lead within this approximation have been 
found in Ref. In the right lead, they are given by 

u{E,x) — uq cosh {a + Ax / d) , (18) 
ip{E,x) = arctan [iJf^"'^ tanh (a -f Ax/c?)] — p. (19) 
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FIG. 2: The density of states N{E, 0) at the device node 

(a) , and the supercurrent spectral density laiE) at = n/2 

(b) , for several values of the ratio L/d, calculated by using 
numerical solution of Eq. 1231 . 

The solution in the left lead are obtained by the change 
of the signs of x, (f), a, and p. The spatial constants in 
Eqs. (|18|l and H19I) can be parameterized as 



7^' 



A = Arccosh 

Uq 



Arccosh^, (20) 

Uq 

uo = coshOo — uocosha — us{E, A), (21) 
Wo = us{E,A), E^Ecosp, A = A cos((/)/2+p),(22) 



where us{E,A) = iJ/V E^ — A^, and expressed via a 
single parameter p, which is to be evaluated from the 
equation following from Eq. H17|l . 



uoA sinp ~ a6a. 



(23) 



The magnitude of the parameter p is controlled by the 
parameter a = d/2L. When a decreases, i.e., the re- 
sistance of the injection lead increases, p turns to zero, 
according to Eq. f^ . and the spectral functions ap- 
proach their values in a closed short SNS junctions. In 
the limit a — 0, DOS has the proximity gap \A^\, where 
A^ — Acos(0/2), and reveals a BCS-like singularity at 
the gap edge, 



N{E,0) ^ Ee{E - \A^\)/ ^ 



(24) 



where Q{x) is the Heaviside step function. The super- 
current spectral density Is{E) spreads over the region 
jA^I < -E < A, and has the singularity at the proximity 
gap edge as wellfSI 



Is{E) 



TT A^eiE~\A4 



2d 



E^ 



-e(A-S). (25) 



In a general case, a 7^ 0, the proximity of the normal 
reservoir leads to finite DOS at all energies, as shown in 
Fig.|2Ia), though it is noticeably suppressed at i? < jA^j 
for L ^ d. The supercurrent spectral density at finite 
d/L extends over the whole subgap region [see Fig.l^fb)], 
while at > A, both Is and Ian turn to zero. Thus, 
in short diffusive junctions with transparent interfaces. 



the supercurrent is carried exclusively by the bound An- 
dreev states confined to the potential well formed by the 
junction. However, this result is only correct to zero ap- 
proximation with respect to the small parameter d/^piS 



B. Opaque interfaces 

The effect of the interface becomes important when the 
interface resistance Rns exceeds the resistance of the nor- 
mal conductor Rn = 2d/u, r = Rns/Rn ^ 1- In par- 
ticular, the magnitude of the Josephson current is deter- 
mined by the Rns rather than hy Rn in the limit r ^ 1. 
At the same time, as we will see below, the suppression 
of the proximity effect is governed by much smaller pa- 
rameter r{d/^o)'^ <C r, and the proximity effect can be 
strong even when r ^ 1. 

A high-resistive interface can be modelled by an effec- 
tive tunnel barrier characterized by its resistance Rns in 
the normal state, which results in the following boundary 
conditions for the Green functions^ a,t x = d. 



aRNsdON = unvs cos (0/2 - i/'w) 
(jRnsI = vnvs sin((^/2 - -0Ar), / = 



usvn, (26a) 
(26b) 



and similar for x ~ —d (the indices N and S refer to the 
normal and superconducting sides of the interface). 

In the limit r 3> 1, the spatial variation of the spectral 
phase is strongly non-homogeneous: the phase drops at 
the barriers and is small in the normal region, ip^ ^ 
along with the spectral current /. The spatial variation of 
the spectral function u is small and can be approximated 
by weakly varying parabolic function. 



u{E, x) 



Uq 



1 + (/3/2) {x/dY 



(27) 



In Eq. H27|l. we neglected the effect of the injection lead, 
assuming its resistance to be larger than RnSi l/a^ r. 
The coefficient /? ^ 1 is to be found from Eq. (|15|) . in 
which the electron-hole dephasing effect has to be taken 
into account because it now becomes comparable with 
the small current-induced depairing. 



(3 = -2i- 



E 



2 2 



(28) 



In Eq. H14|l . we may neglect spatial variations of the inte- 
grand which results in a linear spatial dependence of the 
phase, ip{E,x) « Ix/vq ~ r~^. 

By making use of Eqs. H27|) and H28|) . the boundary 
conditions in Eqs. (|26|l give the equation for the spectral 
functions upj and wjv. 



Us 

UN 



i^E 
UN A 



vs (p 
— cos — 
VN 2 



rip%, 7 = 2r(d/Co)', (29) 



and the expression for the spectral current /, 



VnVs 

1 = sm ■ 



ctRns 2 ' 



(30) 
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FIG. 3: Variations in the DOS (a), and in the supercurrent 
spectral density (b) with the phase difference in SINIS junc- 
tion, calculated from Eqs. and 1321 at 7 = 1: ij> = O.Itt 
(dotted lines), O.Gtt (solid lines), and O.Qvr (dashed lines). 

in which we omitted the small phase ipN from the trigono- 
metric functions. Equation H29I) describes three mech- 
anisms of depairing. The first term in the r.h.s. rep- 
resents the electron-hole dephasing within the normal 
metal. The parameter 7 determines the magnitude of 
the energy gap Eg ~ A/(l -f- 7) in the spectrum of the 
horizontal lead (see below). The second term describes 
suppression of the condensate function vm due to rapid 
change of the spectral phase across the tunnel barrier. 
This effect is similar to the mechanism which produces 
the Andreev bound states in the vicinity of the tunnel 
junction in the ballistic^-^^ as well as diffusive'^^ Joseph- 
son structures. The third term is caused by the supercur- 
rent flow through the normal lead. Neglecting this small 
(^ r~^) term, we obtain the solution of Eq. (|29|l . 



spectral density essentially differs from Eq. H25|) . 

sin(/) A2e(£;- |A0|)e(A-S) 



Is{E) 



2(tRns 



- A2 VA2 - E^ 



(33) 



It is interesting to note that the expression for the equi- 
librium Josephson current obtained from Eq. (|33|l . 



A sin0 
2eRNS 



K{\sm{cj)/2)\), 7«1 (T-O) (34) 



{K is the elliptic integral), can be reproduced with the 
arguments of the scattering theory, similar to the case of 
perfect diffusive SNS junction, by applying the trans- 
missivity distribution for a normal symmetric double- 
barrier structure, -'^ Such a possibility is explained by the 
absence of electron- hole dephasing in this limit. 

The proximity gap is strongly suppressed, Eg <^ A, 
only in the limit of very large barrier strength, 7^1^ 
This is the effect of enhanced electron-hole dephasing, 
similar to the case of a long diffusive SNS junction, where 
the proximity gap is also reduced due to the dephasing 
effect and closes at = tt. This situation is qualitatively 
different from the case of the tunnel junction with a single 
barrier, giving the Josephson current in junctions with 
two strong barriers. 



A sin (/) 
j.. = t:-^: In 



47 



2ei?jvs7 |cos(0/2)p 



7»i (r = o), (35) 



to be much smaller than the result of the tunnel model^ 



A 



VN 



Ve'^-a^' 



A{E, 



l + 7VA2-i;2/A' 



(31) 



According to Eq. , the energy gap Eg (0) in the spec- 
trum of the junction is to be determined by the equation 
Eg = \A{Eg, whose solution can be well approxi- 
mated by a simple relation Eg = |A0|/(1 -I- 7). 

Thus, the regime of strong proximity effect with the 
proximity gap being of the order of A [Fig. 3(a)], per- 
sists in short junction despite of high-resistive interfaces, 
Rns ^ Rn, as long as the condition 7^1 holds. In this 
case, the supercurrent spectral density. 



,Inr/ = -^^Im 



aR 



NS 



Avn{E,^ 



/E2 



, (32) 



extends over all quasiparticle states above the proxim- 
ity gap, including the continuum states above the bulk 
energy gap A, where Is{E) is negative [Fig. 3(b)]. In 
the limit 7 <C 1 (recall that we nevertheless assume here 
the interface resistance to be large, r ^ 1), the dephas- 
ing effect becomes negligibly small and the energy gap 
approaches [A^l, similar to the perfect SNS junction dis- 
cussed above. However, the shape of the supercurrent 



IV. KINETIC EQUATIONS 

In the absence of inelastic collisions, the kinetic equa- 
tions in each lead have the form of conservation laws for 
the spectral currents I±{E), 



D±df± + IsU ± landf^ = I±{E) = const. 



(36) 



At the junction node, the conservation law for the ma- 
trix currents in Eq. (|3J) imposes the boundary condition 



1^ 



(37) 



where the indices L, R, and V refer to the left, right, 
and injection leads, respectively. At the transparent in- 
terfaces, the distribution functions are determined by the 
local-equilibrium population in the reservoirs, 

/+(±rf)=tanh^ iE>A), /_(±d) = 0, (38a) 



tanh ■ 



E 



■eV , E 
± tanh ■ 



eV 



2T 



2T 



,(38b) 



At i? < A, the quasiparticle population in the leads is 
disconnected from the superconducting reservoirs due to 
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complete Andreev reflection, and the quasiparticle den- 
sity function /_|_ is determined by the condition of the 
absence of the net probability current, Ij^ = 0. Due to 
the conservation law in Eq. H37I) . the subgap probability 
current /+ turns to zero within the entire device. 

In the energy region > A, where the currents Is and 
Ian turn to zero, and the diffusion coefficient Z?-|_ turns 
to unity, the kinetic equations have a simple solution. 



no 



i?4 



2R\ 



I^ = - 



2R\ 



(39) 



Here the quantities R+ = d, R^ ~ L, i?_ = d{l/D^'^), 
RX, = L{\/DY.) playthe role of effective resistances of the 
leads for the spectral currents I± , and the angle brackets 
denote spatial averaging along the leads. The currents in 
the left lead are equal by magnitude but flow in opposite 
directions, = —I±, and therefore the currents in the 
injection lead are twice the currents in the left lead. Com- 
bining this result with the relation /+(0) — no = —R+lf , 
following from the kinetic equations, we find that in the 
limit of long injection lead the boundary condition at the 
junction node becomes independent of applied voltage. 



/+(0)==no, L>d, 



(40) 



which implies that the quasiparticles in horizontal leads 
are in equilibrium with the superconducting reservoirs. 

Within the subgap energy region, E < A, the situation 
is more complex due to appearance of the currents Is and 
Ian] the only simplification is due to the zero quasipar- 
ticle current, /+ = 0. By this reason, — const = n+ 
within the entire injection lead, including the junction 
node. Thus, the boundary conditions for the distribu- 
tion functions in the horizontal leads read, 

/+(0) = n+, /_(0)=n_-i?r/r, /-(±rf) = 0. (41) 

Taking advantage of the symmetry of the quantities 
D±{x) = D±{-x), Ian{x) = -Ian{-x), wc Separate 
the even and odd parts of the distribution functions, 
J''''{x) = [f{x) ± /(-a;)]/2, in Eqs. which then be- 
come split in the two independent pairs of kinetic equa- 
tions. One pair that couples and 

D+dfl + Isfl-Iandfl = 0, (42a) 

D^dfl+Isn+Iandfl - (/i^ + /^)/2, (42b) 

has a constant solution, = "-+, fl = 0, consistent 
with the boundary conditions, which yields the relation 
Isn+ = {l^ + I^) /2. As we will see later, Eq. (gHJ, 
the non-equilibrium Josephson current js has the form 
js — {j^ + j^) /2, and taking into account Eq. we 
arrive at the following result jSS, 



Js 



^ r dEIs{E)n+{E). 
e Jo 



(43) 



The second pair of kinetic equations couples the func- 
tions and /I, 

D+dfl + IsZ-Iandfl = 0, (44a) 
D^dfl+Isf^+Iandfl = -lY/2. (44b) 



Since the source term and also the boundary conditions 
to these equations, Eq. 1)41(1 . depend on lY, these func- 
tions determine the dissipative current. The solution to 
Eq. I|44|) in general case must be found numerically. 

At zero temperature, it is possible to further extend 
the analysis. By making use of a step-wise shape of the 
distribution functions n± = 8[±(i? — eV)], we find that 
a trivial solution, /° — 0, ft = 0, lY = 0, satisfies 
Eq. H44|l and all the boundary conditions a,t E > eV. 
Thus the dissipative current vanishes in this energy in- 
terval. On the other hand, ai E < eV, where — 
and n_ = 1, equations 1)44(1 have a non-trivial solution, 
which implies that the dissipative current exists at these 
energies, while the Josephson current is zero, according 
to Eq. ((43(1 . Thus at zero temperature the dissipative 
and non-dissipative currents flow within the separate en- 
ergy regions, which do not overlap: The injection current 
spreads over the energy region < _B < eV, 



f = — dElY{E), 



(45) 



while the supercurrent occupies the region eV < E < A, 



Js 



e{A-eV)— / dEIs{E). 

e JeV 



(46) 



The analysis for the subgap region also applies to the 
case of resistive interfaces, r ^ 1. However, at the ener- 
gies E > A, the supercurrent Is and anomalous current 
Ian are nonzero and give additional contribution to the 
Josephson current in Eq. 146() . 



NONEQUILIBRIUM JOSEPHSON 
CURRENT 



In equilibrium, the Josephson current is given by the 
second term in Eq. (Q, as it was mentioned in Sec. II. Un- 
der non-equilibrium conditions, this connection becomes 
ambiguous and needs reconsideration. The reason is that 
the appearance of the dissipative currents and related 
gradients of the distribution functions in Eq. jS)) will lead 
to spatial variation of the supercurrent term along the 
horizontal lead. To find an appropriate equation for the 
observable non-equilibrium Josephson current, we refer 
to a generic definition of the dc Josephson effect, as a 
current fiow through a junction without any dissipation. 
In our case, the rate of the energy transfer from a voltage 
source to the junction is given by equation. 



1"^ dV 
J Jx j{x)— = j'^V' 



+ j"V 



(47) 



where V'"'^ are the voltage drops at the left/right leads. 
The stationary Josephson effect assumes zero voltage 
drop between the superconducting electrodes, V^+V^' = 
0, thus the non-dissipative current component must sat- 
isfy the equation = = js ■ Combining this equation 
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0.0 0.2 0.4 0.6 0.8 1.0 
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FIG. 4: Josephson current vs voltage at different phases: solid 
lines - L — d, dashed line - L ^ d (a), and Josephson current 
vs phase at several values of the ratio L/d (b). 



with the KirchhofF's rule, we arrive at the following def- 
inition of the Josephson current through the currents in 
the left, right, and injection leads, 



J 



L,B. 



Js±j''/2. 



(48) 



Thus, in order to evaluate observable Josephson current 
in a general case, it is necessary to calculate the injection 
current and the current in one of the horizontal leads, 
and then apply Eq. (|48|1 . In particular case of symmetric 
junction, this procedure leads to Eq. (|^ . 

Persistent current in a SQUID is the most fundamental 
manifestation of the Josephson effect. The Josephson 
current in Eq. H48|l coincides with the circulating current, 
and it can be directly measured by measuring the induced 
flux with an external magnetometer. 

Let us first consider junctions with transparent inter- 
faces, where the non-equilibrium Josephson current is 
given by Eq. (|^ . Since the spectral density Is is pos- 
itive in this case, as it is found from numerical solution 
of Ea. H23|l . and the population of the subgap states is 
depleted with increasing voltage, the injection will sup- 
press the Josephson current and block it completely at 
eV > A [see Fig. ^fa)]; however, the current direction 
cannot be reversed. The Josephson current weakly de- 
pends on the applied voltage and is close to the equi- 
librium value as long as the voltage is smaller than the 
proximity gap value, eV < jA^j. We note that this equi- 
librium value differs from that in closed SNS junctionSf2S> 
it is reduced due to the proximity of the normal reservoir 
and therefore depends on the length of the injection lead, 
as shown in Fig. 0^b). At larger voltage, the Josephson 
current- voltage dependence jsiV) becomes more steep, 
especially for the small phase differences. 

For a long injection lead, Eq. H46|l takes the form. 



b,V) 



dE 



E^-^l 



(49) 




0.0 0.2 0.4 0.6 0.8 1.0 

(])/7l 

a) 



0.0 0.2 0.4 0.6 0.8 1.0 

eViK 

b) 



FIG. 5; Current-phase relations from Eq. 15011 at different 
voltage (a), and critical current vs voltage (b) at L 2> rf. 
For a given voltage, the current-phase relation follows the 
equilibrium law as soon as the proximity gap edge is above 
the energy eV , but it is significantly suppressed when the gap 
edge is below eV . 



explicitly performedjSS 
ttA 



In- 



1 + sin{<i)l2) 



= ,(50) 



2ei?jv /(y) + v//'(V^)-cos2(0/2)' 
/(y) = max[el//A,cos(0/2)], el/ > |A^|. (51) 

This current- voltage dependence is shown in Fig. Efa) 
by a dashed line. In this case, the Josephson current at 
eV < |A^| is constant and equal to the equilibrium value. 

To estimate the efficiency of the Josephson transistor 
let us consider the most steep part of the current-voltage 
characteristic, js{4i: V), at small phase and at large volt- 
age. For example, at (j) = O.Stt, when the equilibrium 
Josephson current approaches about of 0.7 of its critical 
value jc, the switching effect occurs within a small volt- 
age interval 5V ~ O.lA/e. The current gain in this case, 
Sjs/Sj^ ~ 0.7jc/G5V ~ 7{L/d), exceeds unity even for 
comparable lengths of the leads, and it can be further 
enhanced by making the injection lead longer. The up- 
per bound for the gain is imposed by the condition of 
small quasiparticle dwelling time /D, compared to the 
quasiparticle relaxation time r, /D <C t. 

The non-equilibrium Josephson current-phase depen- 
dence for the junctions with transparent interfaces and 
high-resistive injection lead is shown in Fig. |SJa) for dif- 
ferent applied voltages. The kinks on the graphs corre- 
spond to the phase values, at which the applied voltage 
equals the proximity gap, (jioiV) = 2 arccos(ey/A). At 
smaller phases, < (j)< (j)o{V) , the current-phase depen- 
dence has an equilibrium form, while at larger phases 
it is considerably distorted. Correspondingly, the criti- 
cal current jc{V) remains independent of applied volt- 
age until (t)o{V) exceeds the value 0™ = 1.97, at which 
the equilibrium supercurrent approaches its maximum 
value jco = jc(0) = 0.667rA/2ei?. At larger voltage, 
> Acos(0m/2) = 0.55A, the critical current decreases 
and turns to zero at eV > A, as shown in Fig. I^b), 



At zero temperature, the integration in Eq. (|49|) can be 



jc{V) 3s[MV),0] 



Jc{0) 



Jc{0) 



eV A 
1.51— Arccosh— . (52) 
A eV ^ ' 
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FIG. 6: Inversion of the current-phase relation js{(t>) under 
applied voltage in a SINIS junction at 7 = 1. The voltage 
step between the curves in right panel is 10~^A/e. 



In junctions with resistive interfaces, r 3> 1, the cur- 
rent-phase dependence is more interesting, because of the 
possibility of the Josephson current inversion and the 
crossover to the 7r-junction regime. This results from 
the negative contribution of the energies i5 > A to the 
Josephson current which turns to zero before the volt- 
age achieves the gap value, eV < A, when the positive 
and negative parts of Is{E) compensate each other. At 
larger voltage the current becomes negative, as shown in 
Fig-inja). Detailed analysis of the crossover region can be 
made for the junctions with high-resistive injection lead, 
i/d 3> r ^ 1, and at zero temperature. In this case, 
in the horizontal leads, the distribution function /_ is 
small, and the function /+ is approximately constant and 
approaches the equilibrium value nQ{E) — 1 in the su- 
perconducting reservoir [see Eq. ^^]. By these reasons, 
the small dissipative and anomalous components can be 
omitted from the current spectral density, /_ ~ Is, which 
then becomes independent of applied voltage at eV > A. 
This results in the following modification of Eq. 1)46(1 , 



dEIs{E), 



(53) 



min{eV.A) 



where Is{E) is to be found from Eqs. and (|^ . 

As follows from Eq. H53|l. the critical voltage, at which 
the current turns to zero, depends on the phase, and 
therefore the crossover extends over a certain, in fact 
rather small, voltage interval, as shown in Fig. Elb). 
When the voltage approaches the critical region, a new 
current node in the current-phase dependence splits from 
the node at (/< = tt, then it moves towards smaller 
(/>; the process ends when the extra node approaches 
(/) = 0. Such a fine structure of the Josephson current 
inversion has been observed experimentally in long SNS 
junctionsiii At very large interface resistance, 7 1, 
this fine structure becomes irresolvable because of in this 
limit the phase dependence in Is{E, (f>) for relevant ener- 
gies is given by a prefactor sin 0, and therefore the com- 
pensation effect appears simultaneously at all phases^ 



VI. CRITICAL CURRENT IN CURRENT 
BIASED JUNCTION 

In experiment, the current bias setup is often employed 
for investigation of the dc Josephson current. Fig. ^c). 
In equilibrium, the maximum value of the current flowing 
through the junction without creating a voltage drop co- 
incides with the maximum current in the current-phase 
dependence. This is not the case for a non-equilibrium 
junction with current injection: the "critical" current is 
contributed by both the non-equilibrium Josephson cur- 
rent and the injection current. 

Suppose the voltage is applied between the injection 
electrode and left superconducting electrode. Fig. ^c), 
then the external transport current Jt is equal to the 
current in the right lead. In this case, the problem of 
the critical current evaluation is reduced to the analysis 
of the phase dependence of — js — at given volt- 
age. The requirement of zero potential difference between 
the superconducting electrodes is automatically fulfilled 
in our calculation (time- independent phase difference). 
For simplicity, we consider the junction with perfect in- 
terfaces, where the currents are given by Eqs. (|45|) . H46|l . 

The numerical results of such analysis are shown in 
Fig. Ela,c). They are obtained by solving numerically 
Eq. H23|l for the spectral functions and Eq. I|3()|l for the 
distribution functions, which determine the magnitude of 
the injection current. At eV < A, when the supercur- 
rent is allowed to flow through the junction, the current- 
phase relations are similar to that depicted in Fig.[SJa). 
At these voltages, the Josephson current coexists with 
the normal current flowing out of the injection lead. At 
larger voltages, eV > A, the supercurrent is blocked, 
however, the transport current still flows through the 
junction without voltage drop across it, within a certain 
range of the current magnitudes determined by the am- 
plitude of the dependence j^'{(j>). The existence of such 
Josephson-like regime without real Josephson current has 
been first pointed out for a 4-terminal SNS junction with 
opaque interfaces.^ 

To understand this phenomenon, it is important to re- 
member that the injection current in NS interferometers 
is not uniquely determined by the bias voltage, but also 
depends on the superconducting phase. In principle, a 
similar regime with zero voltage drop across the junc- 
tion may appear even for normal reservoirs, at the trans- 
port current jr = {V)/2. This value is unique for the 
given injection voltage, and therefore the corresponding 
dependence jriV) is represented by a straight line. In 
the superconducting junctions, such line broadens to a 
stripe, jT = {V,(f))/2, due to the presence of the free 
parameter cj): The phase adjusts the injection current for 
given injection voltage and transport current to provide 
zero voltage drop across the junction. The width of the 
stripe is determined by the amplitude of the injection cur- 
rent oscillation with the phase. At eV < A, this effect 
is hidden by the presence of the true supercurrent [large 
shaded regions in Fig.[7fb,d)], however, it is fully revealed 
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FIG. 7: Current jV vs phase </> at different injection voltage 
(a,c), and critical transport current jTmm, jr-max vs voltage 
(b,d), at L = d (a,b), and L = 5d (c,d). Shaded regions 
correspond to zero potential difference between the supercon- 
ducting electrodes. In oval-like regions, the Josephson current 
coexists with the normal current; in shaded stripes at eV > A, 
the Josephson current is absent (Josephson-like regime). 

a,t eV > A, where the supercurrent is suppressed [shaded 
stripes in Fig.[3b,d)]. In fact, at large voltage, the width 
of the shaded stripes is determined by the amplitude of 
phase oscillations of the excess injection current. The 
qualitative difference between the phase dependence of 
the excess current (eV > A) and the Josephson current 
{eV < A), is clearly seen in Fig.[7{a,c). It is interesting 
that the "critical current" has different sign for positive 
and negative voltages [the shaded stripes in Fig. [3b, d) 
are differently oriented with respect to the straight line] . 
This is consistent with the fact that the excess current 
changes sign along with the applied voltage. 



VII. INTERFEROMETER EFFECT 

In this Section, we investigate the conductance of the 
injection lead as a function of the bias voltage and su- 
perconducting phase focusing on its properties due to the 
strong proximity effect. At small temperatures, T <^ eV, 
the overall voltage dependence of the differential conduc- 
tance is given by Eq. gHJ, G{V,(j)) = dj^ /dV = alY, 
where the injection current lY is to be calculated by nu- 
merical solution of Eq. 123|) for the spectral functions 
and the kinetic equations H36() . As shown in Fig. [SJ the 
conductance has two peaks, at voltages eV = A and 
eV ~ A^. The first peak is associated with enhanced 
transmissivity of the junction due to the DOS peak at 



FIG. 8: Differential conductance G vs voltage at L = d (a), 
and L = 5d (b) for different phases. 

the bulk gap edge. This peak has a logarithmic singu- 
larity at (j) — 0, and it becomes smeared and decreases 
while (j) departs from zero. The second peak manifests a 
rapid change in the spectral functions in the vicinity of 
the spectrum edge in a short SNS junction (see Fig.[2Jl, 
and it can be interpreted as a resonance transmission due 
to enhanced DOS at the proximity gap edge [A^l. As 
soon as L increases, this resonance becomes more sharp 
because of the singularities in the spectral functions at 
E ~ [A^l become more pronounced, whereas the peak at 
the bulk gap edge, eV = A, decreases and vanishes at 
djL — > 0, as shown in Fig. (S^b). Furthermore, the con- 
ductance exhibits the reentrance effect: G(F) approaches 
the value Gjv = a jidjl -V L) in the normal state both at 
small and large voltages, as it was predicted for NS point 
contactsi^ We notice that the differential conductance 
G{V) deviates from Gn in the short-arm interferometers 
at the characteristic energy of the order of A, in con- 
trast to the long-arm SNS junctions {d 3> £,q), where 
the conductance peak appears at the Thouless energy 
Eti, = hV/{2df « A.^ 

At zero phase difference, </> = 0, the function G{V) can 
be found analytically. In this case. Is = Ian = 0, and 
the function /_ obeys a simple equation D^df- = /_ in 
each lead, with the diffusion coefficients 



7^^'^ = cosh^ 



cosh 



Re0s{E) 



1 + a\x\/d' 
1 + a 

i-y/L 

1 + a 



From Eqs. lll5|l and l(5^ . we obtain 



G{V) = GNvieV) 
rj{E) = z Arctanh z^"'^, z ■ 



(54a) 
(54b) 



(55) 
(56) 



The oscillations of the conductance peak at eV = A 
with the phase can be found from the following argu- 
ments. At this energy, the diffusion coefficient D- turns 
to infinity in the horizontal lead which therefore be- 
comes non-resistive with respect to the normal current, 
i?_(A) = 0. Thus, the differential conductance at eV = 
A is completely determined by the resistance of the injec- 
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FIG. 9: Differential conductance vs phase ai L — d and dif- 
ferent voltage: V = A/e (a), V = 0.7A/e (b, solid line), and 
V = 1.2A/e (b, dashed line). 

tion lead R^{A) = Ltanh6'o/6'o, where Oq = lncot(0/4), 
G^axW = G(A/e,0) ^ Gn In cot ^. (57) 

According to Eq. (|57(l . the peak height approaches 
Giv(l + a) at = TT. At this point, like at (p — 0, the 
spectral densities of the superconducting and anomalous 
current turn to zero, and only the resistances _R_ and 
RY^ are involved in the calculation. Since the condensate 
function becomes completely suppressed in the middle of 
the junction, — 0, the injection lead behaves as a nor- 
mal wire, and therefore the resistance approaches its 
normal value L. Correspondingly, the resistance of each 
horizontal lead coincides with the resistance of a short NS 
junction with the length d, i?_ = dta.nh(Re9s)/ ReOs, 
and therefore G{V) at </> = tt is given by 

In the limit of long injection lead, d/L 0, its resis- 
tance RY. = LtanhRe6'o/Re6'o, 6*0 = Arctanh(A0/£;), 
completely determines the injection current, and the volt- 
age dependence of the differential conductance shown in 
Fig. ISfb) is approximately described by Eq. (|55|l . with 
jA^I substituted for A in the function ri{eV). 

The differential conductance exhibits full-scale 27r- 
periodic oscillations with the phase difference (j) (An- 
dreev interferometer effect). The form of the oscilla- 
tions is qualitatively different for the subgap bias region, 
eV < A, and for eV > A, as shown in Fig. Elb). In 
the latter case, the phase dependence of G has a com- 
paratively simple form, with maxima at = 2?™ and 
minima at (j> = {2n + l)7r. At eV < A, the differential 
conductance approaches minima both at even and odd 
multiples of tt, which reflects the interplay between the 
position and amplitude of the resonance at eV — jA^j. 

VIII. SUMMARY 

We have developed a theory of the non-equilibrium 
Josephson effect in a three-terminal diffusive interferom- 



eter with short SNS junction having the length 2d much 
smaller than the superconducting coherence length ^o- 
We focused on the case of strong proximity effect, when 
the proximity energy gap in the normal region is of the 
order of A. For the junction with transmissive NS inter- 
faces, the density of states N{E) and the supercurrent 
spectral density IsiE) extend over the whole subgap re- 
gion, < i? < A, due to the proximity to a normal 
reservoir, and exhibit a considerable enhancement at the 
energy equal to the proximity gap \A^\ = A|cos(^/2)| 
in the spectrum of a closed SNS junction. The supercur- 
rent spectral density is positive at all relevant energies. 
We demonstrated a possibility of the strong proximity 
effect in a junction with opaque interfaces whose resis- 
tance Rns is much larger than the normal resistance Rn 
of the junction arms. In such case, the suppression of 
the proximity gap, Eg{(f)) w |A0|/(1 -I- 7), is controlled 
by the parameter 7 = {2Rns / RN){d/ (.o)"^ , which could 
be small in short junctions, c? <C ^O; even at large inter- 
face resistance, Rns ^ Rn- In contrast to the case of 
transmissive interfaces, the supercurrent spectral density 
Is{E) is negative above the bulk gap value. 

In three-terminal interferometers, the supercurrent 
generally coexists with the dissipative current flowing 
out from the injection electrode. In such situation, we 
defined the non-equilibrium Josephson current js as the 
non-dissipative component of the current flowing be- 
tween the superconducting electrodes. In symmetric 
junctions, within the subgap energy region, this compo- 
nent coincides with its intuitive representation through 
the supercurrent spectral density and the quasipar- 
ticle population imposed by non-equilibrium injection, 
because the Andreev reflection blocks quasiparticle ex- 
change with equilibrium superconducting reservoirs. In 
junctions with transmissive interfaces the Josephson cur- 
rent becomes completely blocked at eV = A at zero tem- 
perature, while injunctions with high-resistive interfaces, 
the Josephson current undergoes inversion at eV ^ A, 
which spreads over a finite voltage interval. At the ener- 
gies above the bulk energy gap, E > A, the population 
in the junction arms is basically determined by the equi- 
librium population in the superconducting reservoirs. By 
this reason, the Josephson current becomes voltage inde- 
pendent at eV > A. 

We notice that spectroscopy of the supercurrent spec- 
tral density at the subgap energies is possible at zero 
temperature, similar to the tunnel spectroscopy of N{E), 
because the derivative of the Josephson current over ap- 
plied voltage, djs/dV, is proportional to Is{eV). 

The critical current jc of the three-terminal junction, 
defined as the maximum value of the transport current 
jx flowing through the junction without creating a volt- 
age drop, does not coincide with the maximum in js{<P)- 
This is due to the presence of phase dependent injection 
current {4>) which contributes to jc, along with the 
Josephson current, and adjusts its magnitude providing 
zero voltage drop across the junction. At large voltage, 
where the Josephson current is suppressed, the domain 
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of the Josephson-like regime is determined by the am- 
phtude of phase oscillations of the excess current in the 
injection electrode. 

The behavior of the injection current is highly sensitive 
to the quasiparticlc spectrum of the junction and can be 
used to detect the position of the phase-dependent prox- 
imity gap. In particular, the differential conductance of 
the junction with perfect interfaces exhibits sharp peaks 
at the bulk gap value, eV = A, and at the proximity gap, 
eV = lA^I; the latter becomes more pronounced as the 
resistance of the injection lead increases. Furthermore, 
the differential resistance exhibits full-scale oscillations 



with the phase difference; at eV < A, the shape of the 
oscillations becomes rather complex, due to the interplay 
between the position and amplitude of the proximity gap 
resonance. 
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